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Abstract 


NASP 


Presented is a feasibility and error analysis for a hyper- 
sonic flush airdata system on a hypersonic flight experi- 
ment (HYFLITE). HYFLITE heating loads make intrusive 
airdata measurement impractical* Although this analysis is 
specifically for the HYFLITE vehicle and trajectory, the 
problems analyzed are generally applicable to hypersonic 
vehicles. A layout of the flush-port matrix is shown. Sur- 
face pressures are related airdata parameters using a simple 
aerodynamic model. The model is linearized using small 
perturbations and inverted using nonlinear least-squares. 
Effects of various error sources on the overall uncertainty 
are evaluated using an error simulation. Error sources mod- 
eled include boundary-layer/viscous interactions, pneu- 
matic lag, thermal transpiration in the sensor pressure 
tubing, misalignment in the matrix layout, thermal warping 
of the vehicle nose, sampling resolution, and transducer 
error. Using simulated pressure data for input to the estima- 
tion algorithm, effects caused by various error sources are 
analyzed by comparing estimator outputs with the original 
trajectory. To obtain ensemble averages the simulation is 
run repeatedly and output statistics are compiled. Output 
errors resulting from the various error sources are pre- 
sented as a function of Mach number. Final uncertainties 
with all modeled error sources included are presented as a 
function of Mach number. 

Nomenclature 


Acronyms and Initialisms 
FADS flush airdata sensing 

HY-FADS hypersonic flush airdata sensing 

HYFLITE hypersonic flight experiment 

ICBM intercontinental ballistic missile 

MOS model output statistic 
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National Aerospace Plane 

Knudsen’s empirical constant (24.6) 

^ • . 2 
tube cross-section areas, in 

Knudsen’s empirical constant (2.46) 

sonic velocity, ft/sec; also Knudsen’s 
empirical constant (3.15) 

stagnation point pressure coefficient 
pressure coefficient at incidence angle 

tube diameter, in. 

port diameter, in. 

aerodynamic model function 

Mach number function 

Mach number function 

integer roundoff of real value 

empirical coefficient for port damping 
(0.50) 

tube length, in. 
surface length, in. 

effective tubing length, in. (includes port 
effects) 

firee-stream Mach number 
downstream Mach number 
pressure, lb/ft 2 

sampled pressure value, lb/ft 2 

pressure at incidence angle, 0 , lb/ft 2 

stagnation pressure downstream of shock, 
lb/ft 2 

free-stream static pressure, lb/ft 2 
pressure at transducer, lb/ft 2 



from heating, in. 


p o 

pressure at port, lb/ft 2 

9. 

free-stream incompressible dynamic pres- 
sure, lb/ft 2 

*2 

dynamic pressure downstream of shock, 
lb/ft 2 

R 

acoustic impjedance (lbf-sec/ft 3 ) 

R t 

Reynolds number at surface 

s 

Laplace variable (1/sec) 

t 

skin thickness, in. 

T ... 

average temperature in tube, °R (°F) 

T . 

temperature at outer edge of boundary 
layer, °R (°F) 

P surf 

surface temperature, °R (°F) 

T 

Irons 

transducer temperature, °R (°F) 

T 

w 

adiabatic wall temperature, °R (°F) 

V 

transducer volume, in 3 

v +l 

upstream velocity normal to shock wave, 
ft/sec 

V n,2 

downstream velocity normal to shock 
wave, ft/sec 

V 2 

downstream velocity, ft/sec 

V 

OO 

upstream velocity, ft/sec 

X 

longitudinal coordinate, in. 

Y 

lateral coordinate, in. 

Z 

normal coordinate, in. 

a 

angle of attack, deg 

a c 

thermal expansion coefficient, carbon, 
inyin7°R 

a s 

thermal expansion coefficient, steel, 
inyin7°R 

P 

angle of sideslip, deg 

8 

ramp angle, deg 

5 *p 

pressure altitude residual, ft 

V. 

free-stream Mach number residual, ft 

V. 

free-stream pressure residual, ft 

V 

free-stream dynamic pressure residual, ft 

^Ri 

incremental radius of curvature resulting 


Si incremental length change resulting from 

heating, in. 

angle-of-attack residual, deg 
5p angle -of-sides lip residual, deg 

8 0 incidence angle residual, deg 

<|> clock angle, deg 

y ratio of specific heats 

T oblique shock angle, deg 

T| white noise source 

k Knudsen number 

X cone angle, deg 

\i dynamic viscosity, lbf-sec/ft 

6 incidence angle, deg 

p density, lbm/ft 3 

0 standard deviation 

£ damping ratio 

X viscous interaction parameter 

co n natural frequency, rad/sec 

V gradient of vector 

Subscripts 

1 port index 

j iteration index 

oo free-stream parameter 

0 initial condition of parameter 

2 parameter downstream of shock 

Introduction 

The former National Aerospace Plane (NASP) program 
proposed a hypersonic flight experiment (HYFLITE) to 
evaluate the performance of sc ram jet technologies. 1 The 
proposed vehicle has a forebody shape resembling the 
NASP, and was designed to support a scramjet test engine. 
This j)ap)er presents the results of a feasibility study and 
error analysis performed for a hypersonic flush airdata sys- 
tem (HY-FADS) on the HYFLITE vehicle. The study and 
analyses are based on models generally accepted as valid 
within the aerospace community. Some parts of the system 
are not modeled and axe so noted. Although the analysis 
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performed is specific to the HYFLITE vehicle and trajec- 
tory, the methods described and the problems analyzed are 
generally applicable to other hypersonic vehicles. This 
paper addresses feasibility and accuracy issues common to 
most hypersonic vehicles. 

This study presents the physical layout of a candidate 
flush-port matrix, and the effects of various error sources 
on the overall system uncertainty are evaluated using an 
error simulation. In this simulation a realistic set of pres- 
sure data is generated along the trajectory, and various 
error sources are superimposed on the data. Using the sim- 
ulated pressure data (with the various error sources super- 
imposed) for input to the airdata estimation algorithm, 
estimates of the trajectory parameters are generated. The 
effects caused by the various error sources are analyzed by 
comparing the outputs of the estimation algorithm with 
the original trajectory data. To obtain ensemble averages 
of the errors induced by random input error components, 
the simulation is run repeatedly and model output statis- 
tics (MOSs) are generated. Output error magnitudes 
resulting from the various (input) error sources are pre- 
sented as a function of Mach number. The peak Mach 
number analyzed is approximately 13 and the peak alti- 
tude is 120,000 ft This report analyzes only supersonic 
and hypersonic conditions. 

Background 

Because airdata values are used to normalize perfor- 
mance and flight mechanics parameters and to schedule 
control and guidance system gains or to provide direct 
flight control feedback, accurate airdata measurements 
have always been critical for the flight test community. 
Although the requirements for measurement uncertainty 
for hypersonic vehicles have been difficult to firmly estab- 
lish, uncertainties for Mach number of about 0.10 and 
flow incidence angles of about 0.25° have frequently been 
requested. Accuracies on this level are needed for engine 
flowpath and vehicle trajectory optimization 

Historically, airdata measurements were obtained using 
intrusive booms that extended beyond the local flow field 
of the aircraft and measured airmass velocities by directly 
stagnating the flow via a pitot tube at the end. 2 Flow inci- 
dence angles were measured using mechanical vanes 
attached to the probe. Because heating loads induced by 
hypersonic vehicles make it difficult for intrusive mea- 
surement devices to survive, probes are considered 
impractical for hypersonic vehicles. 

To bypass the difficulties with intrusive systems, the 
flush airdata sensing (FADS) system concept, in which 
airdata are inferred from nonintrusive surface pressure 
measurements, was developed in the early 1960’s and 
used to measure hypersonic conditions on the Saturn IVb 
space launch vehicle and X- 15 rocket. The concept was 


further refined during the shuttle reentry airdata sensing 
experiment, which demonstrated that the concept was fea- 
sible for blunt-nosed hypersonic vehicles. Concurrent with 
the shuttle experiment, early aeronautical applications in- 
cluded programs conducted at NASA Dryden Flight 
Research Center on the KC-135 and F-14 vehicles. 4 * 5 Later 
programs at NASA Dryden demonstrated the autonomous 
capability of the system, with fully redundant systems 
being demonstrated on the F-18 High Alpha Research 
Vehicle and autonomous real-time systems demonstrated 
on the F-18 Systems Research Aircraft. 6 - 7 References 8 
and 9 provide an overview of failure detection and fault 
management techniques developed for the real-time FADS 
system. 

The FADS method offers the advantage of inferring 
dynamic pressure, Mach number, and flow incidence 
angles without requiring that the stagnation flow condition 
be created and measured directly. This approach allows the 
heating load at the pressure orifices to be significantly 
diminished for orifices positioned away from the forebody 
stagnation point and makes the FADS method ideal for 
application to hypersonic vehicles. 

Trajectory and Vehicle Description 

HYFLITE achieves high Mach numbers at high 
dynamic pressures using a depressed trajectory. To reach 
this trajectory an intercontinental ballistic missile boosts a 
test article to hypersonic conditions suitable for perform- 
ing scram jet research. 1 This trajectory accelerates the vehi- 
cle to approximately Mach 13 at 120,000 ft with peak 
dynamic pressures of 2800 lb/ft 2 . Figure 1 shows the pro- 
posed flight envelope. 
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Fig. 1. HYFLITE vehicle flight envelope. 

Figure 2 pictures the vehicle to be analyzed. The pro- 
posed configuration resembles the NASP forebody and 
was designed to support a scramjet test engine. The vehi- 
cle has a shovel-like nose with 5° wedge angles at the 
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Front view 
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Fig. 2. Proposed pressure FADS port locations for HYFLITE vehicle. 


leading edge. The top surface ramp expands to a 6° half 
angle approximately halfway aft, and the lower surface 
has an isentropic compression ramp for the scramjet fix- 
ture. For this study the isentropic ramp was approximated 
by a 6° wedge. The longitudinal leading-edge radius of the 
nose is 0.1 in., and the lateral leading-edge radius is 
approximately 20 in. The approximate vehicle length is 
33 ft. Although the leading edge — which was chosen to 
give good flow characteristics at the scramjet inlet for 
single-stage-to-orbit performance — is quite sharp, this 


analysis assumes that the leading edge is blunt enough to 
allow sufficient heating relief for the nose to withstand the 
resulting surface temperatures without active cooling. 

Airdata System Port Layout 

The port matrix layout was determined as a compromise 
between the need to accurately model the flow characteris- 
tics near the pressure ports, and the need to give good sen- 
sitivity of the pressure measurements relative to the airdata 
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parameters. Regions where high-temperature viscous 
interactions cause boundary-layer thickening were 
avoided. Only locations along the leading edge (where the 
boundary-layer thickness is small) and locations farther 
aft on the ramp surfaces (where the viscous interactions 
diminish) were considered. 


The angles of attack and sideslip, a and p, are related to 
the port coordinate angles 6 by 

cos(0.) = cos (a) cos (p) cos (\) + 
sin (P) sin (0.) sin (X f ) + sin (a) cos (P) cos (<)>,.) sin (X ) 


Results presented in references 6, 7, 8, and 9 show that 
more ports always lead to more accurate measurements, 
but that diminishing returns are reached when the number 
of ports exceeds nine. Below this number of ports, airdata 
measurement errors can rise significantly. Thus, a total of 
nine ports was used: five along the leading edge and four 
distributed on the centerline of the two-dimensional ramp 
surfaces. The ports along the leading edge provide the pri- 
mary information source for measuring Mach number 
(Af), angle of sideslip (P), and static pressure (pressure 
altitude, Hp). The ports on the upper and lower ramp sur- 
faces provide the primary information source for angle-of- 
attack (a) measurement. 

The port locations are described using a clock and cone 
angle coordinate system that is an analytical mapping 
from the actual surface to a hemisphere. The cone angle, 
X , is the total angle that the normal to the surface makes 
with respect to vehicle axis of symmetry. The clock angle, 
<(> , is the clockwise angle looking aft around the axis of 
symmetry starting at the bottom of the fuselage. As men- 
tioned earlier, to keep the locations of the leading-edge 
ports far enough forward on the vehicle so that boundary- 
layer effects could be ignored, ports were placed near the 
nose with cone angles of 0°, ±45°, and ±80°. The ramp 
ports were placed farther aft on the vehicle where the 
effects of viscous interactions diminish. The first set of 
ports were located approximately 4-ft aft of the nose, and 
the second set 8 ft aft. Figure 2 depicts the HYFLITE 
vehicle and approximate port locations. Table 1 lists the 
coordinate angles of the pressure ports. 


Table 1. Pressure port coordinates. 


Port no. 

h deg 

<j>, deg 

Location 

0 

0 

0 

Stagnation 

1 

85 

0 

Lower ramp 

2 

85 

180 

Upper ramp 

3 

84 

0 

Lower ramp (aft) 

4 

84 

180 

Upper ramp (aft) 

5 

45 

-90 

Leading edge 

6 

45 

90 

Leading edge 

7 

80 

-90 

Leading edge (aft) 

8 

80 

90 

Leading edge (aft) 


where 0 ( is the total incidence angle relative to the local 
surface. 


Pressure Modeling 

The key to this analysis is to develop a pressure model 
that relates the airdata to the surface pressures in an accu- 
rate way, but is still simple enough to be inverted and 
implemented as part of an airdata estimation algorithm. 
For accuracy equilibrium high-temperature gas effects 
must be modeled in the analysis. Because Mach number 
peaks out in the Mach 12 region, however, previous 
results 10 indicate that nonequilibrium high-temperature 
gas effects are not significant, and they were ignored in 
this analysis. Because pressure is a mechanical quantity 
and depends primarily upon the mechanical aspects of the 
flow, influences caused by high-temperature gas chemistry 
are always secondary. Viscous boundary-layer interactions 
are significant, however, errors resulting from the interac- 
tions were modeled and their effects introduced in the 
error analysis section. 

Near the leading edge the stagnation properties (C P 
were computed using the high-temperature equilibrium 
normal shock tables and equations of state for air, ll> 12 and 
the pressure distribution was described using modified 
Newtonian flow. 10 






COS 2 (0) = 


Cp cos 2 (6) 

r max 


(la) 


On the ramp surfaces, the pressure distribution was com- 
puted using the equilibrium normal shock tables and the 
shock wave angle— turning angle relationship for high-tem- 
perature equilibrium flow. 10 Here, the oblique shock angle, 
T, is an implicit function of the free-stream velocity and 
the surface wedge (turning) angle, 6 


y n2 

tan(r-S) = r 7 ^tan(D (lb) 

«. 1 

where V n j is the upstream velocity normal to the shock 
wave, and * s *e downstream velocity normal to the 
shock wave. The ratio 
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V .2 V'jSin (H 
V+i " ^sin(D 

is determined from the equilibrium normal shock tables. 
For a given turning angle, equation (lb) is iterated until 
the velocity ratio is consistent with the normal shock 
tables. Once the shock angle and downstream velocity 
have been determined, the downstream pressure coeffi- 
cient is determined from the equations of continuity and 
momentum, 13 


P 2 V 2 ” P\ V v 

P 2 + p 2 [ V 2 sinr] 2 = P x + p, [VjSinr] 2 


to give 


CP wedge ~ ^ S * n ^ 



(1C) 


The specific form of the nonlinear function F[...] 
depends on whether the port is located on the leading edge 
or on the ramp surfaces. The parameter t\ i is a noise vec- 
tor and can be used to represent the various modeling-, 
system-, and measurement-error sources. Taken together, 
the matrix of ports forms an overdetermined nonlinear 
model in that the number of observations (9) is larger than 
the number of states to be estimated (4). 

The system is nonlinear in a and p, and Mach number is 
a parameter of the system, which is in turn related to the 
airdata parameters. As a result of this nonlinearity and 
interdependence, the system of equations must be solved 
iteratively to estimate the airdata states. This iteration is 
implemented by linearizing about a starting value for each 
port, and the perturbations between the measured (gener- 
ated) data and the model predictions are evaluated. At a 
given data frame for each pressure port, i = 1, ... N, 

SP j i + l = M 2 . l\A-, ♦?"-)] = 


The clock angle is always either 180° for the upper sur- 
face ramp or 0° for the lower surface ramp; thus, the inci- 
dence angle is related to the wedge half-angle by 

Upper Surface: 

cos(0 f ) = [cos (a) sin (8) - sin (a) cos (8) ] cos (P) 


Lower Surface: 

cos (0;) = [cos (a) sin (8) + sin (a) cos (8) ] cos (P) 




8pJ + 1 

8a y+1 

8p ; + 1 


+ ... 


Defining 


(2b) 


Equations (la), (lb), and (lc), along with the high- 
temperature equilibrium normal shock tables for air, can 
be used to simulate pressure data along the vehicle flight 
trajectory. This system of equations was also used as the 
model that relates the measured pressure to the airdata 
states in the airdata estimation algorithm. 

Airdata Estimation Algorithm 

This section develops the regression algorithm that was 
used to extract the airdata estimates from the pressure 
data. For this analysis the airdata state vector is described 
in terms of four parameters: dynamic pressure behind the 
shock (< 7 2 ), angle of attack (a) , angle of sideslip ( P), 
and free-stream static pressure (PJ. Using these four 
basic airdata parameters, other airdata quantities of inter- 
est may be directly calculated. For a given pressure obser- 
vation the form of the model is 


+ 

1 


4 


* J+i 
5 <?2 

rj" 

= 

P j 

+ 

8P i +1 

y+i 

a 


a* 


8a >+1 

P > + 1 


_P ; . 


6p y+1 _ 


the updated state vector is solved using recursive least 
squares 


LH 


*2' 


r -i 

^2 

pj +1 

= 

P 1 

OP 

+ [(vV)(VF)] ( V?F ) 

8P/ + 1 

a' + 1 


a' 


8 V +I . 

P y+1 


. _ 


(2c) 


P^a.) = F 



<h 




, <)>,, X p y, M" 


a 






+ 11 ; (2a) 


At the end of each iteration cycle the downstream Mach 
number, Af 2 , is determined from the ratio of q 2 //*<». Poo is 
used to determine the altitude, and the temperature is 
determined using a standard atmosphere. With these val- 
ues for pressure, temperature, and Af 2 , free-stream Mach 
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number is determined using a table lookup of the equilib- 
rium normal shock tables for air. For each data frame the 
iteration cycle is repeated until algorithm convergence is 
reached — typically from 2 to 5 cycles. For this algorithm 
free-stream Mach number is used as the convergence cri- 
terion, and changes in the free-stream Mach number of 
less than 0.001 are considered to indicate a converged 
cycle. At the beginning of each new data frame, the sys- 
tem of equations is linearized about the result of the previ- 
ous data frame. The Choleski factorization techniques 
used to perform the regression are standard imple- 
mentations 14 and will not be discussed here. 


The subscript w indicates quantities at the wall, and the 
subscript e represents quantities at the outer edge of the 
boundary layer. If pressure is constant (vertically) 
throughout the boundary layer and viscosity is propor- 
tional to the square root of the temperature, 11 then 



Assuming an adiabatic wall (with a recovery factor of Rc ), 


Error Analysis: Development 
of Error Models 

In this section error sources for the pressure measure- 
ments are considered, and error propagation models, both 
random and systematic, are developed. These error models 
are used to contaminate the pressure data. The corrupted 
pressure data are subsequently used as input to the airdata 
estimation algorithm. The simulation is designed to allow 
errors to be superimposed individually or in groups. Impo- 
sition of the errors individually allows the magnitudes of 
the various effects to be quantified individually. 


AT 




1/4 


(3c) 


A value of X < 3 indicates a weak boundary -layer inter- 
action, and a value of x > 3 indicates a strong interaction. 
For a flat plate the increment in pressure resulting from a 
weak interaction is approximately 10 


8P 1 

VISCOUS I 

p 

_ inviscid 


1 +0.31x + 0.05x 2 (4a) 


j weak 


This error analysis assumes that the flow models 
described earlier are accurate. Except those made for 
viscous boundary-layer interaction, no attempt is made to 
characterize the uncertainty of the surface pressure 
models. The error models developed address only the 
local measurement error sources. Vehicle dependent 
effects, such as bending about the center of gravity and 
position error, 2 are not considered. For a real vehicle, these 
effects must be calibrated. 

Viscous Interaction Error Model 


and the increment in pressure resulting from a strong inter- 
action is approximately 

[ 5/> v--l » 0.5 14x- 0.241 (4b) 

L ^ inviscid J strong 

If the wall is nonadiabatic such as for an actively cooled 
surface, then the viscous interaction errors will be less 
than those predicted by equations (3) and (4). Thus, the 
viscous interaction analysis is conservative. 


As mentioned previously, equilibrium high-temperature 
gas effects are included in the estimation model, and non- 
equilibrium high-temperature effects are considered negli- 
gible in this analysis. Because the shock is detached near 
the nose, shock-boundary-layer interactions can be 
ignored. At high free-stream velocities, however, viscous 
interactions between the boundary layer and the external 
flow can significantly increase the boundary-layer thick- 
ness, which can in turn significantly affect the local pres- 
sure distribution. Based on the analysis of reference 10, 
the induced pressure error (caused by viscous interactions) 
is directly proportional to the governing similarity param- 
eter, x. where 


M Ip Li 

- OO r w r* w _ 

^ /n . A II 


A/ 


MM J*e 


KT w \l w 


RT 


(3a) 


Pneumatic Lag Error Model 

Because of the hostility of the surface environment, sur- 
face pressures must be sensed by transducers located 
remotely within the aircraft body. The measurement trans- 
ducers are connected to the surface via lengths of pneu- 
matic tubing. This tubing induces spectral distortion 
caused by frictional losses and wave resonance. Based on 
the analysis described in reference 7, the pneumatic atten- 
uation effects are modeled by a time-varying, second- 
order system of the form 


Pl(S) = 1 

p 0 ( s ) J_/ + 2 JL j+ i 

® 2 a. 


(5a) 
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where P L is the pressure measured at the transducer and 
Pq is the surface pressure input. The natural frequency of 
the system is described by 


LV 


" 2 ~ 

" l 

PO 

[l+ iLR)2 ] 

4 _ 

[l + (^) 2 1 

1 

1 

[ 2V J 

12 (p 0 c) 2 


6(p 0 c) 2 

J 


and the damping ratio is given by 

. R(0 n \ L 2 LV 

2 P 0 [ 2 C 2 A e c\ 


where L is the tubing length, A c is the tube cross-sec- 
tional area, c is the local sonic velocity, p 0 is the local 
density (within the tube), and V is the entrapped volume of 
the transducer. Assuming laminar flow in the tube the 
acoustic impedance of the system is defined by 

D 

For each pressure sensor £ and o) n are evaluated at each 
time point based on an average of the surface conditions 
and the conditions at the transducer. 

The effects of the surface pressure ports are modeled 
using the analyses of references 16 and 17. Here the effect 
of the port is to extend the effective length of the tube by 
the amount 


this opposing flow is that a static pressure gradient is 
established, with the cold region of the tube reading lower 
than the hot region. This effect, referred to as transpira- 
tion or molecular creep, is primarily a steady-state 
phenomenon. 

For high ambient pressure levels at which continuum 
flow exists, the molecular creep effect is not considered 
important and can be ignored. At low-pressure levels, 
however, this effect becomes important. For low-pressure 
conditions the pressure gradient induced by longitudinal 
temperature gradients is primarily a function of Knudsen’s 
number, k. If the characteristic length of the system is 
taken to be the tube diameter, then for air 


K- 




\lJT 

D P 


The transpiration effect is modeled by numerical inte- 
gration of Knudsen’s molecular creep equation 18 


dP(x) _ 
dT(x) ~ 

l P(x) r K(p,T,L>) (K(p, 7 M)) +a) I 2 
2F(x) [k (P,T,D) 2 + ( a + b)K(P,T,D ) +bc_ 

(6) 

where the constants a, b, and c have been empirically 
determined to be 


a b c. 


Knudsen’s constants 24.6 2.46 3.15 


g D i rn 

* L .fr-fr+-£r R 'y < 3b > 

where y, P.M , p., and c are evaluated based on local con- 
ditions, and D t is the port diameter. Entrapped vol- 
umes found wnere the tube joins the surface port are 
ignored. For this analysis the effective increment in tubing 
length is small, about 1 to 5 percent of the total length. 

Thermal Transpiration Error Model 

The transpiration-induced pressure nonequilibrium, 
which results from large temperature gradients and small 
tubing diameters, is primarily a molecular flow phenome- 
non. 18 Along an unequally heated gas boundary, kinetic 
theory predicts that the gas adjacent to the boundary wall 
will creep from the colder region to the hotter region. For 
static conditions in which there is no net mass flow within 
the tube, to balance this creep the gas in the center of the 
tube tends toward the colder end of the tube. The result of 


and D is the tube diameter. For a prescribed value of D 
equation (6) may be treated as an initial value problem, 
which is integrated over a prescribed temperature range 
from the transducer to the surface. It can then be imple- 
mented as a part of the overall simulation scheme. 

Port Misalignment Error Models 

This analysis models the effects of port alignment 
errors. TVvo types of errors are modeled: (1) an initial (ran- 
dom) error in the port alignment and (2) systematic defor- 
mation of the structure resulting from heating loads. 

Initial Port Misalignment Error Model 

The clock and cone angle bias alignment errors were 
introduced into the simulation at the beginning of each 
Monte-Carlo time history run using a normalized random- 
number generator. Once the values of the alignment errors 
have been set for a particular data run, they are held con- 
stant throughout the data run. 
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Structural Deformation Error Model 

Modeling of the thermal expansion of the highly curved 
leading edge is extremely difficult, and only the deforma- 
tion of the ramp surfaces is addressed in this model. The 
analysis is for a flat plate with the primary deformation 
mechanism assumed to be differential expansion of the 
skin layers. Referring to figure 3, the skin is assumed to be 
carbon-carbon on the exterior bonded to a stainless steel 
substrate on the interior. The thermal expansion coeffi- 
cients are a c (carbon-carbon) and (^(stainless steel). 
Assuming identical thickness for the two layers, then 
based on the analysis of reference 16, the radius of curva- 
ture resulting from differential thermal expansion of the 
two layers is 


stainless steel), T i is the temperature at the surface, and 
T re j is the reference temperature at which the skins were 
bonded. Referring to figure 3, the differential incidence 
angle is given by 


5 


J Si 



(7b) 


where 8 Si is the length expansion to the surface. From ref- 
erence 17 the expansion per unit length is given by 


*Si 


a 




(7c) 


R = — 

Ri ~ 3(a c -a t ) (T r T ref ) 


(7a) 


where 5^. is the induced radius of curvature at the ith 
port, t is the thickness of the skin (carbon-carbon and 


and the resulting deformation model is then given by 

5., - CM) 



Fig. 3. Thermal expansion deformation model for HYFLITE skin. 
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For this analysis the per-length parameter, L i , is the dis- 
tance over which the heating has occurred and for this 
analysis is defined as the distance from the tip of the vehi- 
cle to the port 

Resolution and Transduction Error Models 

System quantization levels of 12 and 16 bits are ana- 
lyzed. To simulate this quantization the generated pressure 
data are sampled and truncated after n bits (2 n counts) of 
resolution. 


exceeded. Furthermore, it is assumed that known bias 
errors are eliminated using preflight tare calibrations; 
unknown bias errors or variable but systematic transduc- 
tion errors are lumped with the random errors for this anal- 
ysis. The random transducer errors are modeled as 
unbiased Gaussian white noise sequences and are gener- 
ated using a normalized random-number generator. The 
errors are allowed to vary randomly for each data frame 
and across the collection of pressure ports. 

Results and Discussion 


j- p. -\ full Scale Reading „™ e results of the error analysis will now be presented. 

"/ = Int b ii'c — i n ^ x 2 x a Figures 4 through 10 show sample residuals (from one of 

u cae e in 2 the Monte-Carlo data runs) for each of the error sources. 

(8) For the final uncertainty analysis of the system, all error 
sources are superimposed simultaneously. The model out- 
This analysis assumes that the measurement transducers puts from 100 data runs are averaged to give ensemble s la- 

are located in an environmentally controlled area, and that tistics. Table 2 lists root-mean-square (1 -ct rms) error 

the manufacturer’s environmental specifications are not statistics. 
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Fig. 4. Mach-number residual caused by high- 
temperature viscous-boundary-layer interactions. 



(a) Mach-number residual. 

Fig. 5. Residuals caused by pressure-tubing lags, reso- 
nances, and attenuations. 




(b) Angle-of-attack residual. 
Fig. 5. Concluded. 



(a) Mach-number residual. 

Fig. 6. Residuals caused by thermal transpiration of 
pressure tubes. 
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(b) Angle-of-attack residual. 
Fig. 6. Concluded. 



(a) Mach-number residual. 



(b) Angle-of-attack residual. 

Fig. 7. Residuals caused by initial port misalignment 
error. 



(a) Mach-number residual. 



(b) Angle-of-attack residual. 

Fig. 8. Residuals caused by beating-induced structural 
deformation. 



(a) Mach-number residuals. 


Fig. 9. Residuals caused by transducer error. 
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(b) Angle-of-attack residuals. 
Fig. 9. Concluded. 



(a) Mach-number residuals. 

Fig. 10. Residuals caused by measurement system reso- 
lution error. 



(b) Angle-of-attack residuals. 
Fig. 10. Concluded. 


Errors Caused by Viscous Interactions 

In this analysis the data were simulated using the analy- 
sis of equations (3) and (4) to superimpose viscous interac- 
tion pressure errors onto the baseline data Because the 
analysis of equations (3) and (4) is invalid near the leading 
edge, errors introduced by boundary-layer interactions 
were modeled only for the ramp ports and the 80° pressure 
ports. The estimation algorithm was run, and the resulting 
trajectory was compared with the original trajectory. The 
Mach-number residuals between the estimated airdata and 
the actual trajectory were computed, and figure 4 shows 
these results. The effect on Mach number is significant 
only beyond Mach 8, with a peak error of 0.17 (at Mach 
12.5). Because the upper and lower ramp pressures are 
affected equally by the viscous interactions, little influence 
is experienced by the angle-of-attack measurements. 

Errors Caused by Pneumatic Lag and Attenuation 

In this section the errors induced by the pneumatic 
arrangement of the airdata system are analyzed. Because 
the damping ratio and the time lag of the pressure mea- 
surement sensors vary inversely with the mean density 
level in the tube, the responses of the ports near the stagna- 
tion point are much less damped. In addition the lags for 
these ports are significantly shorter than those for the ports 
located away from the stagnation point where the pressure 
and density levels are lower. Because the HY-FADS algo- 
rithm is time independent, this mixture of time delays and 
damping manifests itself as a data distortion that becomes 
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Table 2. 1-a rms absolute and 1 -<j rms percentage error estimates for the HY- 
FADS system. 



V 

5-, 

lb/ft 2 

Up' 

ft 

««. 

deg 

6p, 

deg 

% 

6--. 

q/q 

% 

°Hp/Hp> 

% 

1.0 

0.0 

1.0 

10 

0.2 

0.1 

0.3 

0.1 

0.1 

2.0 

0.0 

1.5 

10 

0.2 

0.2 

0.2 

0.1 

0.0 

3.0 

0.0 

0.3 

10 

0.1 

0.2 

0.1 

0.0 

0.0 

4.0 

0.0 

0.8 

25 

0.1 

02 

0.1 

0.0 

0.1 

5.0 

0.0 

4.0 

100 

0.2 

0.3 

0.2 

0.2 

0.2 

6.0 

0.0 

5.0 

250 

0.2 

0.3 

0.4 

0.6 

0.2 

7.0 

0.0 

2.0 

350 

0.2 

0.3 

0.6 

0.3 

0.3 

8.0 

0.1 

1.0 

450 

0.2 

0.4 

1.0 

0.1 

0.4 

9.0 

0.1 

3.0 

500 

0.2 

0.3 

1.1 

0.4 

0.5 

10.0 

0.1 

3.0 

600 

0.2 

0.3 

1.2 

0.3 

0.5 

11.0 

0.2 

2.0 

750 

0.2 

0.3 

1.6 

0.2 

0.7 

12.0 

0.3 

1.0 

1100 

0.2 

0.3 

2.3 

0.1 

0.8 

12.5 

0.4 

4.0 

1350 

0.2 

0.3 

3.3 

0.4 

1.0 


most pronounced at high Mach numbers, at which the 
spread of damping and lags is greatest Thus, it is gener- 
ally a good idea to match the lags and damping of the var- 
ious sensors as closely as possible over the range of test 
conditions. A practical system with minimal pneumatic 
effects (except an inherent time-latency) can generally be 
achieved by varying the lengths and diameters of the con- 
nective tubing. Table 3 lists the configuration used for this 
simulation. 


Table 3. Pneumatic layout of airdata system. 


Port 

no. 

X, 

deg 

<t>. 

deg 

Tube 

diameter, 

in. 

Tube 

length, 

in. 

Port 

diameter, 

in. 

0 

0 

0 

0.020 

60 

0.010 

1 

85 

0 

0.060 

24 

0.020 

2 

85 

180 

0.060 

24 

0.020 

3 

84 

0 

0.060 

24 

0.020 

4 

84 

180 

0.060 

24 

0.020 

5 

45 

-90 

0.020 

60 

0.010 

6 

45 

90 

0.020 

60 

0.010 

7 

80 

-90 

0.060 

24 

0.020 

8 

80 

90 

0.060 

24 

0.020 


For this analysis it is assumed that the transducers are 
located in a thermally controlled environment and that the 
transducer temperature is maintained at 100 °F for the 
duration of the flight The volume entrapped by the trans- 
ducer is assumed to be minimal, and a value of 0.05 in 3 
was used. To eliminate problems with surface heating, the 
external port size is constrained to be 0.01 in. for all lead- 
ing-edge ports and 0.02 in. for the ramp surface ports. The 
model parameters are computed using the average of the 
surface and transducer temperatures. Temperatures at the 
external ports are computed using the equilibrium oblique 
shock equations for the ramp ports, and the equilibrium 
normal shock tables for the leading-edge ports. The gas 
density within the pressure tubes was evaluated using the 
input pressure level and the average between the surface 
temperature (computed from equilibrium analysis) and the 
assumed transducer temperature (100 °F). Figure 5(a) 
shows the resulting Mach number residuals, and the angle- 
of-attack residuals are shown in figure 5(b). 

Errors Caused by Thermal Transpiration 

The thermal transpiration or molecular creep error anal- 
ysis was performed using the configuration analyzed in the 
previous section. Figure 6 presents Mach-number and 
angle-of-attack residuals. Again the effect on Mach num- 
ber is significant only beyond Mach 6. It is at these speeds 
where the conditions rarefied enough to induce molecular 
effects. The effects on angle of attack are insignificant. 
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Errors Caused by Port Misalignment 

This section presents the results of the port misalign- 
ment error analyses. Errors resulting from the initial port 
misalignment will be presented first The errors resulting 
from systematic deformation of the structure caused by 
heating loads will be presented next 

Initial Port Misalignment Error 

For this analysis the three initial misalignment error 
bounds were considered — 0.05°, 0.10°, and 0.25°. The 
squares of these bounds defined the variances for the 
random-number generator. Figure 7 shows the results of a 
single data run. 

Thermal Expansion Port Misalignment Error 

Figure 8 presents the effects of the thermal expansion 
on port misalignment As mentioned earlier, the modeling 
of the thermal expansion of the leading edge is beyond the 
scope of this analysis, and only the deformation of the 
ramp surfaces was addressed. The reference temperature 
was assumed to be 70 °F and the skin thickness 0.5 in. The 
surface temperature was computed using the equilibrium 
oblique shock equations. A value of 4.4 xl0^/°R was 
used for a c , and a value of 7.562 x 10^/°R was used for 
a s . Computations were based on the distance from the 
nose — which is unconstrained and can move freely — to 
the pressure port. 

Errors Caused by Transduction and Resolution 

The transduction uncertainty and system resolution are 
important and must be carefully considered when devel- 
oping a HY-FADS system. Depending on the needs of the 
measurement users, varying degrees of raw measurement 
quality and sampling resolution is required. For this analy- 
sis it was assumed that absolute pressure transducers were 
used. The analysis assumed a 40 lb/in 2 (5760 Ib/ft 2 ) full- 
scale reading for the three transducers connected to the 
ports nearest the stagnation point. A 15 lb/in 2 (2160 lb/ft 2 ) 
full-scale reading edge was assumed for the rest of the 
transducers. Three transducer error levels were analyzed; 
these represent 

1. Conventional sensors with 1-a standard error at ±0.5 
percent of full scale 

2. High-quality sensors with standard error at 1-a 
±0.05 percent of full scale 

3. State-of-the-art sensors with standard error of 1-a 
±0.01 percent of full scale 

Figure 9 shows the effects of the transduction error. The 
Mach-number error for the conventional sensors exceeds 
0.96 at Mach 12; for the high-quality sensors the 
Mach-number error exceeds 0.36 at Mach 12. The 


state-of-the-art sensors give a Mach-number error less 
than 0.06 at Mach 12. 

Once the data were corrupted with measurement noise, 
they were sampled assuming 12-bit and 16-bit resolutions. 
Figure 10 presents the effects of the resolution errors. The 
Mach-number error induced by the 12-bit resolution is 
nearly 0.50 at Mach 12. The 16-bit resolution gives an 
error that is less than 0.05 at high Mach number. The effect 
on angle of attack is small for both resolution levels. 

Collected Error Analysis: Monte-Carlo Simulation 

Accuracy estimates for a state-of-the-art system were 
assessed. In this system all of the error models were incor- 
porated. The system was assumed to have 16-bit resolu- 
tion, and the transducers were assumed to have a noise 
level of less than 0.01 percent of full scale. The port layout 
was according to table 1 and the initial port layout uncer- 
tainty was assumed to have a 1-a error of 0.05°. The 
sensor pneumatics used the configuration of table 2. To 
obtain ensemble averages of the errors induced by random 
input error components, the simulation was run repeatedly 
and model output statistics (MOSs) were generated. As 
such the simulation can be referred to as a Monte-Carlo 
simulation. 

The 1-a rms absolute error MOSs were listed for Mach 
number, dynamic pressure, altitude, angles of attack and 
sideslip, and as a function of free- stream Mach number in 
table 3. Presented also were the 1-a rms percentage (of 
reading) error MOSs for Mach number, dynamic pressure, 
and altitude. As mentioned earlier, these results are contin- 
gent on the assumption that the aerodynamic models used 
are accurate. 

Even at very high altitudes and Mach numbers at which 
the thermal conditions become the most extreme, the 1-a 
measurement uncertainty still was less than 0.1 for at 
Mach 8, and less than 0.4 at Mach 12. The dynamic pres- 
sure measurement uncertainty was less than 5 lb/ft 2 for all 
of the Mach numbers analyzed. 

Concluding Remarks 

This paper addressed feasibility and uncertainty issues 
associated with applying a flush airdata sensing (FADS) 
system to the hypersonic flight experiment (HYFLITE) 
vehicle, which is a sharp-nose configuration of the type 
required for single-stage-to-orbit performance. Heating 
loads to be experienced by the HYFLITE make the use of 
intrusive airdata measurement devices impractical. 
Although the analysis in this paper was performed specifi- 
cally for the HYFLITE vehicle, the methods described and 
the problems analyzed are considered generally applicable 
to hypersonic vehicles. 
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This study presented the physical layout of a candidate 
flush-port matrix, and the effects of various error sources 
on the overall system uncertainty were evaluated using an 
error simulation. In this simulation a realistic set of 
pressure data was generated along the trajectory and vari- 
ous error sources were superimposed on the data. The 
pressure data contaminated with the various error models 
were then used as input to an airdata estimation algorithm 
and effects resulting from the various error sources were 
analyzed by comparing the outputs of the estimation algo- 
rithm with the original trajectory data. Ensemble model 
output statistics were generated, and output error magni- 
tudes resulting from the various (input) error sources were 
presented as a function of Mach number. 

The error sources modeled include errors introduced by 
boundary-layer thickening caused by high-temperature 
viscous interactions, pneumatic lag induced by the pneu- 
matic layout of the measurement sensors, and thermal 
transpiration resulting from high-temperature molecular 
creep in the pressure sensor pressure tubing. Other error 
sources modeled are misalignment errors in the initial 
pressure matrix layout, port misalignment errors resulting 
from thermal warping of the vehicle noseshape, sampling 
resolution error, and transducer sensing error. Nonequilib- 
rium high-temperature gas effects and vehicle-dependent 
effects, such as bending about the center of gravity and 
position error, were not considered. 

Simulation uncertainty estimates for a state-of-the-art 
system were presented. Even at very high altitudes and 
Mach numbers, at which the measurement conditions 
become the most extreme, the 1 -a measurement uncer- 
tainty still was less than 0.1 for Mach numbers up to 8, 
and less than 0.04 at Mach 12. The measurement uncer- 
tainty of dynamic pressure was less than 5 lb/ft 2 for all 
Mach numbers. 

While the results presented are contingent on the 
assumption that the aerodynamic models used are accu- 
rate, a large range of potential error sources was modeled, 
and the results presented are believed to represent the 
types of accuracies that can be achieved by a hypersonic 
flush airdata sensing (HY-FADS) system. This report indi- 
cates that the HY-FADS measurement system is not only 
feasible for the HYFLITE vehicle, but that without 
extraordinary measures it is possible to obtain sufficiently 
accurate airdata results. The pressure measurements can 
be taken with very small diameter (0.01 in.) tubes to avoid 
problems with surface heating. A flush airdata system 
should be considered for future hypersonic test vehicles. 
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